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I.  Introduction 


The  analysis  of  electromagnetic  pulse  (EMP)  internal 
coupling  to  critical  electronic  components  and  subsystems 
is  complicated  by  the  presence  of  many  seemingly  random 
parameters,  such  as  the  relative  posit ipns  of  bunched  cables 
near  points  of  entry  (POE)  and  the  random  positions  of 
conductors  in  N-wire  lines.  These  random  parameters  pre- 
clude the  strictly  deterministic  solution  for  EMP  induced 
excitations  at  particular  load  points.  One  can,  of  course, 
choose  to  analyze  a single  deterministic  "average  model" 
of  the  system  in  the  hope  that  the  excitations  obtained 
will  indicate  expected  excitations  on  any  of  several  randomly 
different  actual  systems.  If  the  random  parameters  strongly 
affect  the  coupling  to  certain  critical  system  points  the 
actual  excitations  may  differ  vastly  from  the  deterministic 
predictions. 

For  the  case  of  strong  random  effects  in  the  system 
a statistical  analysis  should  be  performed  in  orde'^  to  obtain 
a valid  range  of  expected  excitations. 

The  purpose  of  this  report  is  to  give  a brief  dis- 
cussion concerning  a possible  method  for  the  statistical 
analysis  of  the  load  excitations  on  an  unshielded  N-wire 
random  cable  illuminated  by  an  incident  monochromatic  field. 
The  technique  utilizes  the  concepts  of  time-harmonic  electro- 
magnetic field  reciprocity  and  statistical  representation  of 
an  ensemble  by  a subset.  Although  the  discussion  here  is 
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restricted  to  a limited  class  of  structures  (e.g.  unshielded, 
unbranched  N-wire  cables) , the  method  should  be  extendable 
to  shielded  and  branched  cables  as  well.  In  addition,  it 
may  be  possible  to  conduct  the  analysis  directly  in  the  time- 
domain  using  Prof . Welch' s time-domain  reciprocity  theorem, 
Ref.  (1). 


II.  Probabilistic  System  Description 


Consider  an  ensemble  of  possible  N-wire  cable  config- 
urations, whv^re;  one  such  typical  structure  is  shown  in  Fig.  1 
for  the  case  of  N=3.  The  wires  are  assumed  to  run  in  the 
general  direction  of  the  z-axis.  The  random  relative  posi- 
tion of  the  K-th  wire  in  a constant-z  cross  section  is  given 
by  the  vector  Rj^(z),  as  shown  in  Fig.  1. 

The  probabilistic  system  description  can  be  formally 
represented  as  a stochastic  random  process  where  the  random 
vector  is  actually  a function  of  two  parameters,  s and  z, 
where  s represents  points  in  an  abstract  sample  space,  S, 
having  a 1-1  relationship  to  the  ensemble  of  possible  cable 
configurations.  To  each  subset  of  sample  points  (e.g.  each 
subset  of  possible  cable  configurations) , there  is  assigned 
a positive  probability  via  the  mapping  into  the  real  line 
segment  (0,1)  as  illustrated  in  Fig.  2. 

The  joint  probability  distribution  function  of  the 
positions  of  the  wires  is,  in  general,  quite  complicated. 
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Figure  1.  Random  Wire  Locations 
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Figure  2.  Probability  Mapping 
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It  depends  upon  all  values  of  z since  the  position  of  each 
wire  is  affected  by  the  positions  of  all  wires  everywhere 
on  the  line.  The  usual  approximation  that  is  made  to  sim- 
plify the  statistical  description  of  such  a system  is  to 
assume  that  the  random  vectors  (in  this  case  the  wire 
positions)  are  independently  distributed.  This  allows  one 
to  define  the  joint  probability  distribution  function 
(p.d.f)  as  the  product  of  the  p.d.f.s.  of  each  wire. 

The  assumption  of  independence  may  be  quite  erroneous  for 
the  case  of  a tightly  bunched  cable  system  and  will  not 

be  employed  in  the  statistical  models  to  be  presented. 

A common  method  for  describing  an  ensemble  of  param- 
eter dependent  random  vectors  (e.g.  a stochastic  process) 
is  through  a random  coefficient  basis  function  expansion, 

as  outlined  in  Ref.  (2) ., 

\(z|s)  = <»>j(2)[aj^(s)x  + b.^(s)y]  (D 

j=l 

A 

Where  ^^j(z)j*  is  a complete  basis  function  set,  x and 
y are  the  usual  unit  vectors,  and  |ajj^(s)  , 

coefficient  functions  of  the  sample  point,  s,  in  the  param- 
eter space  related  to  the  ensemble  of  possible  cables. 

The  probability  distribution  functions  of  the  random 
coefficients  will  depend  upon  the  basis  set  used  and  the 
type  of  cable  configuration,  such  as  straight  wire  parallel 
cable,  helical  twisted  wire  cable,  or  twisted  pair  cable. 
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etc.  Given  the  type  of  cable,  it  should  be  possible  to 
deduce  the  form  of  the  p.d.f.  for  the  coefficients. 

A second  method  for  describing  the  random  variation 
of  the  cable  configuration  is  via  a random  walk  process. 

Ref.  (4).  In  this  type  of  description  we  approximate  the 
continuous  possible  wire  locations  by  a discr':tized  model 
as  shown  in  Figs.  3 and  4 having  Q possible  cable  positions 
for  the  Q cables  to  occupy  on  a mutually  exclusive  basis. 

At  each  periodic  cross  section,  where  jumps  are 
allowed  there  are  Q1  possil  cable  arrangements.  At 
the  Z2  cross  secticn  the  probability  of  each  cable  arrange- 
ment is  given  in  terras  of  the  Qi>«Q!  transition  matrix,  P(Z), 
where  = Prob  ^the  ith  configuration  being  realized  at 
Z|  given  the  jth  configuration  at  Z-AZ 

III . Use  of  the  Reciprocity  Theorem 

To  obtain  the  induced  currents  in  the  loads  of  a 
transmission  line  that  is  illuminated  by  some  specified 
spatially  distributed  and  oriented  driving  field  one  can 
solve  either  the  original  boundary  value  problem  or  the 
reciprocal  problem, as  in  Ref.  (3). 

As  an  example  of  the  use  of  the  reciprocity 
theorem,  consider  an  arbitrary  N-wire  line  shown  in  Fig. 

5 being  illuminated  by  the  discrete  electric  dipole 
source,  5 (r-r^)  , where  |Jg|  »1. 
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Figure  3.  Stepwise  Changes  in  Cable  Model 
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ifjure  -1.  Hunched  Cable  Configuration  in  a Cross  Section 
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To  obtain  the  induced  voltage  across  a load  impe- 
dance, say  for  example  Z^,  we  place  a reciprocal  unit 
current  source,  J2>  across  and  compute  the  reciprocal 
source  generated  electric  field  at  the  original  source 
point,  r^.  Then,  using  the  Lorentz  reciprocity  theorem 
for  a region  containing  all  ssource  distributions, 

dv  (2) 

where  Mj^=M2=^»  we  obtain  the  induced  load  voltage  via 


f Ej-Jj+Ej-Mj  dv  = J 

v"  v" 


V 


L 


(3) 


For  the  case  of  a unit  magnitude  discrete  magnetic 
dipole  illuminating  source,  we  can  determine  the 
induced  current  in  Z_  by  placing  a reciprocal  discrete 
voltage  source,  as  shown  in  Fig.  6,  in  series  with  Z^ 
and  determining  the  generated  at  r^.  The  load  current 
induced  in  Z^^  by  is  given  by 

I,  = hI-mT 

The  use  of  the  reciprocity  theorem  to  obtain  induced 
load  impedance  currents  and  voltages  requires  the  ability 
to  compute  the  fields  radiated  by  the  arbitrary  N-wire  line 
when  the  line  is  excited  at  the  load  position  by  a discrete 
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J,  = J,  S(f-r,) 
ILLUMINATING  SOURCE 


mZL  (t)RECIPROCAL 


I'iqure  5.  Arbitrary  N-Wire  Line  with  Reciprocal  Excitations 


Fiquro  6.  Reciprocal  Voltage  Source 
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current  or  voltage  source.  This  may  be  quite  a formidable 
task,  expecially  if  the  N-wire  is  located  close  to  other 
obstacles  such  as  arbitrarily  oriented  cables,  bulkheads, 
etc.,  whose  scattering  properties  must  be  considered. 

IV.  Minimum  Subset  Representation 

To  obtain  useful  knowledge  about  the  expected  induced 
load  currents  and  their  variances  we  must  consider  in  detail 
the  induced  load  currents  in  a finite  subset  of  the  continuous 
parameter  ensemble  sample  space  ^seS^.  That  is,  using  the 
reciprocity  theorem,  we  must  compute  the  radiated  fields 
generated  by  localized  sources  at  the  load  positions.  This 
can  usually  be  done  in  two  steps  by  assuming  weak  coupling 
to  nearby  scatterers  (except  massive  ground  plane  structures) . 
For  the  weak  coupling  case  we  can  first  compute  the  currents 
induced  in  the  N-wire  line  by  the  localized  sources,  neglect- 
ing the  effects  of  the  scatterers.  We  then  consider  the 
field  source  to  be  the  N-wire  structure  radiating  in  the 
presence  of  the  scatterers  and  hopefully,  using  approxima- 
tions, we  can  obtain  at  least  an  upper  bound  on  the  radiated 
field  at  the  illuminating  source  point.  We  will  only  con- 
sider here  the  first  step  in  this  procedure. 

In  addition  to  considering  the  structures  as  being 
random,  the  incident  field  may  also  have  to  be  treated  as 
random.  Some  reasonable  bounds  on  the  maximum  expected 
EMP  induced  excitations  are  the  desired  quantities  of  any 
internal  interaction  analysis,  j-f  the  maximum  drive  incident 
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EMP  field,  {i.e.  spatial  distribution  and  field  orienta- 
tion) , is  not  known  a priori  for  a particular  load  point 
under  consideration  then  the  statistical  expectations  and 
variances  of  load  responses  will  have  to  be  averaged  over 
all  possible  incident  fields  as  well  as  possible  structure 
variations. 

We  must  compute  the  induced  currents  on  a represen- 
tative set  of  the  ensemble  of  N-wire  lines  (possibly  over 
a ground  plane)  generated  by  the  localized  reciprocal  sources. 

A major  question  of  great  importance  is:  "what  constitutes 

a representative  subset  of  the  given  ensemble  of  N-wire 
lines  being  considered?"  The  detailed  answer  to  that  ques- 
tion will  have  to  be  found  at  a later  date  but  will,  of 
course,  depend  upon  the  type  of  cable  configuration.  At 
this  time  we  will  only  consider  some  simplifications 
incurred  by  assuming  certain  subsets  to  be  valid  representatives 
of  the  total  en.emble  for  statistical  averaging  purposes. 

The  simplest  possible  subset  of  random  N-wire  cables, 
as  shoitfn  in  Fig.  1,  is  the  set  of  uniform  straight  parallel 
wires  with  random  distances  between  wires.  The  transmission 
line  currents  on  the  cable  can  easily  be  computed  in  terms 
of  the  random  vector  distances  between  cable  wires.  These 
random  variable  currents  will  then  set  up  random  variable 
fields  which  can  be  obtained  from  a straightforward  computa- 
tion, Ref.  (5).  Using  multiple  transformations  of  the  known 
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probability  densities  of  the  v;ire  locations  a density  on 
the  radiated  fields  can  be  obtained  and  the  average  and 
standard  deviations  of  the  fields  can  be  obtained,  at  least 
in  principle.  The  density  of  the  physical  locations  of 
the  wires  in  this  simple  model  will  usually  be  uniform  for 
all  positions  within  a given  cylindrical  boundary.  The 
advantage  of  using  this  subset  of  random  N-wire  cables  is 
that  in  most  cases  simple  uniform  transmission  line  theory 
can  be  used  to  obtain  the  reciprocal  source  induced  line 
currents.  However,  the  expectation  and  variance  of  the 
load  excitations  obtained  using  this  simple  model  may  not 
correctly  represent  the  corresponding  statistics  for  the 
entire  random  cable  ensemble. 

Another  possible  ensemble  subset  is  the  random 
coefficient,  expansion  given  in  (1)  . The  series  will  have 
to  be  trucated  and  the  random  coefficient  density  functions 
will  need  to  be  deduced.  If  the  maximum  longitudinal 
oscillation  of  the  random  line  is  slow  enough  the  non- 
uniform  transmission  line  theory  can  be  applied  to  obtain 
the  line  currents  in  terms  of  the  random  coefficients.  A 
direct  solution  for  line  currents  can  also  be  obtained 
using  a coupled  integral  equation  approach  where  the  random 
wire  positions  appear  in  the  kernels. 

The  random  walk  description  of  the  cable  wire  posi- 
tions can  also  be  used  to  generate  another  ensemble  subset. 
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I In  this  case  a series  of  piecewise  uniform  cables  can  be 

generated  using  a step  by  step  numerical  Monte  Carlo  random 
walk  procedure.  Each  of  the  realized  cables  is  then  analyzed 
using  multiple  section  transmission  line  theory. 

In  using  the  reciprocity  theorem  the  random  wire 
positions  on  the  cable  will  not  only  play  a part  in  deter- 
mining the  reciprocal  source  generated  line  currents  but 
will  also  appear  as  independent  parameters  in  the  solution 
for  the  random  radiated  field. 

V.  Summary 

A brief  and  somewhat  informal  discussion  has  been 
presented  concerning  the  statistical  analysis  of  load 
impedance  excitations  induced  on  a random  N-wire  cable 
! by  an  incident  time-harmonic  field.  The  N-wire  unshielded 

and  unbranched  random  cable  is  modeled  by  a stochastic 
random  process  which  can  be  approximately  described  using 
a variety  of  common  techniques.  The  use  of  the  reciprocity 
theorem  is  described  for  incident  fields  generated  by 
discrete  electric  or  magnetic  dipole  sources  and  its 
applicability  is  directly  extendable  to  distributed  sources. 
A qualitative  discussion  is  given  with  regards  to  the 
subset  representation  of  a statistical  ensemble  sample 
space.  The  motivation  for  this  concerns  the  possible 
reduction  in  statistical  modeling  complexity  by  considering 
only  certain  geometrically  simple  subsets  of  a complex 
random  system. 

( 

I 

i 12 


The  statistical  modeling  approach  applied  to  large 
and  complex  random  cable  structures  is  an  extremely  impor- 
tant topic  with  respect  to  EMP  internal  interaction  analysis. 
Considerable  work  needs  to  be  done  to  develop  practically 
applicable  methods  for  use  on  real  cable  systems.  The 
emphasis  here  has  been  on  a nonrigorous  presentation  of 
ideas  which  it  is  hoped  will  serve  as  a precursor  for  a 
more  detailed  study  to  be  initiated  in  the  near  future. 
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